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Abstract 

We have simulated (3 + 1)— dimensional finite temperature Z 2 gauge theory by using 
Metropolis algorithm. We aimed to observe the deconfinement phase transitions by 
using geometric methods. In order to do so we have proposed two different methods 
which can be applied to 3— dimensional effective spin model consisting of Polyakov 
loop variables. The first method is based on studies of cluster structures of each 
configuration. For each temperature, configurations are obtained from a set of bond 
probability (P) values. At a certain probability, percolating clusters start to emerge. 
Unless the probability value coincides with the Coniglio-Klein probability value, 
the fluctuations are less than the actual fluctuations at the critical point. In this 
method the task is to identify the probability value which yields the highest peak 
in the diverging quantities on finite lattices. The second method uses the scaling 
function based on the surface renormalization which is of geometric origin. Since 



this function is a scaling function, the measurements done on different-size lattices 
yield the same value at the critical point, apart from the correction to scaling terms. 
The linearization of the scaling function around the critical point yields the critical 
point and the critical exponents. 
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1 Introduction 

Lattice studies of quantum chromodynamics (QCD) show that there exist two dif- 
ferent phase transitions at two extreme limits. The first one is the infinitely heavy 
quark mass limit. This limit is relatively easy to study since QCD with static quarks 
can be studied as pure gauge theory. In this limit, gauge theories undergo a phase 
transition from a low temperature phase in which quarks and gluons are confined, to 
a high temperature phase in which colored quarks and gluons exist. This transition 
is named as the deconfinement phase transition and can be studied in terms of its 
order parameter, Polyakov loop (L) pQ, which is similar to the spin variables in fer- 
romagnetic systems. The deconfinement phase transition of SU(N) and Zjy gauge 
theories are well-understood in terms of Svetitsky and Yaffe (S-Y) conjecture j2]. 
According to S-Y conjecture, (d + 1)— dimensional finite temperature SU(N) and 
Zjy lattice gauge theories are in the same universality class as d— dimensional 
spin system. The problem appears when one wants to study gauge theories with 
finite-mass quarks. In the finite temperature lattice gauge theories, the quark fields 
couple to Polyakov loops similarly to the coupling of external magnetic field to spin 
variables in the ferromagnetic systems [3]. Hence the Polyakov loop is no longer 
the order parameter for the deconfinement phase transition. There have been some 
recent attempts to device new observables for studying the deconfinement phase 
transition in the finite quark mass limit [3J The importance of the definition 
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of a new order parameter is that it may enable one to study the relation between 
the deconfinement and chiral phase transitions which is seen in the second ex- 
treme limit where the quarks are massless. It is known that quarks have finite and 
relatively low mass. One task is to understand whether the deconfinement phase 
transition continues into finite- mass region or not. 

The relation between ferromagnetic phase transition and the percolation theory has 
been established some years ago by Coniglio and Klein jHllZj. One of the proposed 
solutions to the above mentioned problem is based on the studies of the percolating 
clusters existing in the system around the phase transition region. This approach 
is used to introduce a geometric method of following the deconfinement phase tran- 
sition in lattice gauge theories [HI El HOI HH EES 1131 El- The main difficulty in 
applying Coniglio-Klein approach to lattice gauge theories is to obtain a finite effec- 
tive spin Hamiltonian from lattice gauge theory action. In finite temperature lattice 
gauge theories Polyakov loop is the unique gauge invariant quantity which is sen- 
sitive to symmetry breaking. From (d + 1)— dimensional finite temperature lattice 
gauge theory, one can obtain a d— dimensional effective spin system whose dynamic 
variables are Polyakov loops [21 E3- in this work we aim to apply two different 
techniques to the effective spin system obtained from the finite temperature gauge 
theory. The success of these techniques are previously tested on spin systems for 
first- and second-order phase transitions [T71 ITS} FEU I2TH l2T| 12^]. 
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The first approach which we will apply to the (3 + 1)— dimensional finite temperature 
Ising gauge theory is the method which is based on determining the Coniglio-Klein- 
like probability value which produces the largest fluctuations in the 3— dimensional 
effective spin system consisting of the Polyakov loop variables [17]. The second 
approach that we want to implement on lattice gauge theories uses the scaling 
function based on the surface renormalization [T8J HH1 EE EQ 122] ■ Both methods 
can easily be implemented to follow the deconfinement phase transition into finite 
mass quark region. 

The paper is organized as follows. In Section |2] we explain how to implement the 
approaches to finite temperature lattice gauge theories. In Section |3] we have pre- 
sented our simulation results and finally, Section |U is devoted to the conclusions and 
further possible applications of the methods. 

2 Model and the Method 

The finite temperature on the lattice is defined by introducing finite and periodic 
temporal direction, while in all d spatial directions, the lattice is infinite. For com- 
puter simulations this can be implemented by taking the spatial extent (N a ) of the 
lattice much larger than the temporal extent (N T ). In this work in order to dis- 
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cuss the operators of geometric origin to study deconfmement phase transition, we 
have considered the simplest possible gauge theory, namely Z 2 lattice gauge theory 
During our simulations the temporal extent of the lattice is fixed as N T = 2. The 
Hamiltonian of the Ising gauge theory is 

H = £ U^i)Uv(i + + v)Uj(j) (1) 

where U^{i) represents gauge field which takes only two values, U^{{) = ±1, in this 
model. The Polyakov loop variables are given by 

L = ( 2 ) 
For each spatial lattice site, one can assign a Polyakov loop variable. This 3- 
dimensional reduced system can be considered as a spin system. For finite tem- 
perature Z2 gauge theory this 3— dimensional spin system has dynamic variables 
L = ±1. 

In this work we propose to use two different methods with geometric origin. The 
first method is based on the introduction of a set of probability values which will 
determine the Coniglio-Klein-like clusters for each inverse temperature (3 J7j. The 
probability value, which coincides with the actual Coniglio-Klein probability at the 
transition temperature of the effective spin system, can be determined by studying 
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the highest fluctuations in the system. The largest fluctuations are the indication 
of the phase transition ^7j. Hence, one expects that the thermodynamic quantities 
obtained by using this probability value obey the scaling rules giving the critical 
exponents of the system. 

For completeness we will give a brief description of this method. In a ferromagnetic 
system each configuration contains clusters of like spins. These clusters percolate 
before the phase transition point and it is known that these geometric clusters are 
not sufficient to establish the connection between the percolation of these geometric 
structures and thermal phase transitions. By using Fortuin-Kasteleyn clusters 
Coniglio and Klein [7] showed that the partition function of the Ising Model can be 
written in purely geometric terms. In this description of the system one uses a bond 
probability value to define the number of spins which belong to the same cluster. For 
the Ising Model this bond probability for like spins is given as P = 1 — e ~ 2J / kT , This 
probability reduces the existing geometric clusters of like spins in each configuration 
to smaller clusters which are called Coniglio-Klein clusters. It is shown that in all 
dimensions, the thermal phase transitions can be studied in terms of percolating 
Coniglio-Klein clusters. In order to obtain Coniglio-Klein clusters, in other words, 
to find the correct bond breaking probability, one needs to know the Hamiltonian 
of the system. In the finite temperature lattice gauge theory case this Hamiltonian 
is the Hamiltonian of the effective spin system. Assuming that the probability that 
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breaks the existing geometric structures are unknown. In this case, one can choose 
a set of trial probabilities between and 1 and decide which one does correspond 
to the correct probability value. Simulations are performed for a set of chosen 
temperature values near the transition temperature. For each temperature, a set of 
predetermined trial probability values may be used to break the existing geometrical 
connections between like spins at each configuration. In fact this range is not very 
wide since the percolation threshold at each dimension is known for any given lattice 
geometry Moreover a simple test run can limit the range to obtain required 
accuracy. If the original configuration has a spanning cluster of like spins, a subset 
of probability values may result in spanning clusters of Coniglio-Klein type. As one 
uses the bond-breaking probability P, the tuning of this parameter may point out the 
connection between the dynamics of the system and this artificial parameter. As this 
parameter approaches to the value given by Coniglio and Klein at a given inverse 
temperature (3, these clusters coincide with the dynamically generated Coniglio- 
Klein {C-K) clusters and P value can be identified by the largest fluctuations for 
each f3 . Particularly, the value of the parameter P which corresponds to f3 c (L) 
results in the most widely fluctuating cluster sizes. The fluctuations in the system 
can be followed by considering any fluctuating quantity. In this work we have 
chosen the fluctuations of the largest cluster in each configuration as the fluctuating 
quantity. 
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For a certain probability P the average of the largest cluster is defined as, 



1 Nc 1 

<c max > P = w j2 v (cr x (p)) (3) 

c i 

where N c is the number of configurations, C^ ax (P) is the largest cluster obtained 
by using the probability P in i th configuration and V is the volume of the system. 
The fluctuations (MCF(P)) of the largest clusters is given as 

MCF(P) =< C 2 max > P -< C max > 2 P . (4) 

These observables have the characteristic scaling behavior. < C max > P scales like 
magnetization and MCF(P) scales like susceptibility; 

<C max>P = Nj»- d C(tN?) (5) 

and 

MCF(P) = Nl YH ~ d MCF{tN^ T ) (6) 

where N a is the linear size and d is the dimension of the system, t is the reduced tem- 
perature, Yjj and Yt are the magnetic and thermal critical exponents respectively. 
In this method for each simulation temperature, the fluctuations of the largest clus- 
ter will be studied for a predetermined set of probability values. The P value which 
results in the highest fluctuations for a given (3 value is taken as the value of the 
Coniglio-Klein probability value. Particularly the highest peak for the simulation 
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temperatures indicates the transition temperature. 



The second method uses the scaling function F based on the surface renormalization. 
The function F is studied in detail for the Ising Model fHl El 120] and g-state Potts 
Model [2HI22]- I n order to calculate this function, one considers the direction of 
the majority of spins of two parallel surfaces which are L/2 distance away from 
each other. In the calculations, a counter is increased by one only if two opposing 
surfaces have majority of the spins in the same direction. For the Ising Model, the 
formulation of this scaling function is simple, since the direction of the majority of 
the spins can be calculated as the sign of the sum of the spins in a surface [2D] . 
Hence for the a spin system consisting of Ising spins, F can be written in the form 



where Si is the sum of the spins in the i surface. Scaling function F obeys the 
scaling relation 



If F is expanded in Taylor series near the critical point (t — > 0) in the form 



F =< sign [Si] sign [S i+L/2 ] > 



(7) 




(8) 




(9) 



S 



the function F can be written in the linearized form. By using Eq.(jOJ) one can obtain 
the thermal critical index as well as the phase transition point of the effective spin 
system. 

The next section is devoted to our simulation results. 

3 Results and Discussions 

We have simulated (3 + 1)— dimensional Ising gauge theory by using Metropolis 
algorithm. In our simulations we have discarded 2 x 10 4 iterations for thermalization. 
The errors are calculated by using Jack-Knife analysis on the measurements of 10 
to 20 different runs of 10 5 to 10 6 iterations after thermalization. Calculations are 
repeated for 40 (3 values and for 15 P values. In our simulations the time extent of 
the lattice is kept fixed, N T = 2. We have used four different size spatial lattices, 
N a = 10, 12, 16 and 20, but we have presented the results for the largest three lattice 
sizes. 

We have presented our results in two subsections. In the first subsection (Method I) 
the results of the computations of fluctuations in the system with respect to varying 
probability values are presented. In the second subsection (Method II) the results 
obtained from scaling function based on the surface renormalization are given. 
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3.1 Method I 

This method is based on studies of cluster structure of each configuration by repeat- 
edly restructuring the configurations obtained from a set of bond probability (P) 
values. For each lattice size, 15 P values are used in order to obtain Coniglio-Klein- 
like clusters. In this way the criticality of the system is searched by studying cluster 
distribution at each inverse temperature (5. For each (5 value at a certain probabil- 
ity, percolating clusters start to emerge. Nevertheless, unless the probability value 
coincides with the C — K probability value, the fluctuations are less than the actual 
fluctuations at the critical point. The task is to identify the probability value which 
yields the highest peak in the diverging quantities on finite lattices. 

In Figure la we have presented the fluctuations (MCF) of maximum size clusters on 
lattice N a = 20. In this figure we have plotted the maximum fluctuations obtained 
at different temperatures. Here one can see that the fluctuations increase as the 
temperature approaches to a certain value (expected finite-size critical value for the 
effective spin system) and reaches to its peak. The temperature corresponding to the 
peak value is considered as the transition temperature. In Figure lb three highest 
fluctuation peaks are plotted for closer inspection. 

Figure 2a shows the highest fluctuation peaks for spatial sizes N a = 12, 16 and 20. 
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In figure 2b we have presented the scaling of cluster fluctuations. The position of 
the peak yields the phase transition point /3 c (N a ) for finite lattice size N a and the 
peaks of the fluctuations scales like susceptibility with the lattice size in the form 
N% YH ~ d . The f3 c (N a ) is obtained by minimizing the distance between scaled cluster 
fluctuation curves. The thermal and magnetic critical exponents obtained after 
this minimization as Y H = 2.498 ± 0.002, Y T = 1.587 ± 0.002 and /3 c (N a = 12) = 
0.4256 ± 0.0003, (3 c (N a = 16) = 0.4254 ± 0.0002, (3 c (N a = 20) = 0.4253 ± 0.0002. 
The difference between the inverse transition temperature /3 t (oo) for the infinite 
lattice and !3 t (N a ) obtained for different size lattices is given by 

fc{N a ) - fr(oo) ~ N; y t . (10) 

By using Eq. ffTUj) and calculated values of c (N a ) and Yt we have calculated the 
infinite lattice critical point as /3 c (oo) = 0.4251 ± 0.0001. 

3.2 Method II 

The scaling function F given by Eq.((7|) is also based on geometric structures on 
the lattice. In Figure 3a we have plotted the scaling function obtained from 
3— dimensional effective spin model for three different-size lattices in the full (3 range. 
In Figure 3b we have presented the linearized scaling function in Eq. Q around the 
critical point for all three different-size lattices. Since F is a scaling function, differ- 
ent size lattices yield the same value at the transition point apart from the correction 
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to scaling terms. From the pairwise crossing and the slopes we have obtained phase 
transition point j3 c and the thermal critical exponent Yt using Eq.Q. Calculated 
critical exponents and the transition point values are given in Table 1. These two 
methods yield similar results. Particularly, the location of the transition point by 
using two different methods are in good agreement. 

4 Conclusions 

In this work we have aimed to demonstrate two methods to follow the deconfine- 
ment phase transition into finite mass region of lattice quantum chromo dynamics. 
Starting from the gauge theory action one can obtain an effective spin Hamiltonian 
(at high temperature limit) |15j . By using this effective Hamiltonian, in principle, 
one can calculate various quantities to search the phase transition. Nevertheless 
using full effective action is both cumbersome and subject to approximations. In 
the literature this approach is pursued only for the pure gauge theory case and also 
spin systems with external magnetic field. Avoiding the complication of an effective 
Hamiltonian one can still obtain accurate information on the behavior of the system 
through the deconfinement phase transition. None of the two methods described in 
the present work requires the knowledge of the interactions between the Polyakov 
loops. Moreover, since only spin-like dynamic variables are under consideration, the 
interactions between matter fields and the gauge fields do not effect the conclusions 
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drawn by using these two methods. For (3 + l)-dimensional finite temperature Ising 
gauge theory, we have obtained critical exponents Yu and Yr from Coniglio-Klein 
like clusters and Y T by using scaling function F. The results obtained by using both 
methods are accurate and satisfactory for conclusive evidence for the universality 
class of the phase transition. 
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Table Captions 

Table 1. The critical exponents and the transition point values obtained by using 
F function. 
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Figure Captions 

Figure 1. The fluctuations MCF of maximum size clusters on lattice N a = 20. a) 
for all PmaxiP), b) for three P ma x((3) values with highest peaks. 

Figure 2. a) Highest fluctuation peaks for spatial sizes N a = 12, 16 and 20. b) The 
scaling of cluster fluctuations given in a). 

Figure 3. a) The scaling function F for three different size lattices N a = 
12, 16 and 20. b) Linearized scaling function F around the critical point 
for the same lattices. 
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N — N 


Y T 


P 


12 - 16 


1.568 ±0.002 


0.4251 ±0.0002 


16-20 


1.666 ±0.006 


0.4251 ±0.0002 



Table 1: 
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